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, Abstract. We present high angular resolution interferometric observations [ ] of the dust 

continuum at 2.7 and 1.3 mm, and of the HC3N (J=12-^l 1) and C^*0 (J=2^1) emission 



around L1157-mm, a Class object that drives a spectacular molecular outflow. The mil- 
limeter dust emission is clearly resolved into two components, a flattened compact source 
of -^450x250 AU at 1.3 mm, and mass '^O.l Af©, plus an extended envelope of ~3000 AU 
at 1.3 mm, and mass ~1.1 Mq. The millimeter spectral index varies throughout the region, 
with the lower value found toward the compact protostar, possibly indicating grain growth 
in the denser regions. A strong interaction between the molecular outflow and the close pro- 
tostellar environment is taking place and affects the structure of the innermost parts of the 
envelope. This is shown by the spatial coincidence between the molecular outflow and the 
dust (1.3 mm continuum) and HC3N emission: both tracers show structures associated to 
the edges of the outflow lobes. Basically, the global picture sketched for the Class object 
LI 157-mm by Gueth et al. (1997) is supported. We find possible evidence of infall, but we 
do not detect any velocity gradient indicative of a rotating circumstellar disk. 

Key words, individual: LI 157 — stars: formation — stars: circumstellar matter — ISM: dust 
— ISM: molecules radio lines: molecular 



* Based on observations carried out with the IRAM Plateau de Bure Interferometer. IRAM is supported 
by INSU/CNRS (France), MPG (Germany) and IGN (Spain). 
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1. Introduction 

Class low-mass young stellar objects are deeply embedded in circumstellar dust and gas ma- 
terial, and are found associated with very energetic molecular bipolar outflows. Theory outlines 
a scenario where a central object is surrounded by an infalling envelope that contains most of 
the mass (e.g. Larson 1969; Adams, Lada, & Shu 1987). This infalhng material is accreted onto 
the central protostar furmeled through a circumstellar disk that grows as the system evolves (e.g. 
Cassen & Moosman 1981; Terebey, Shu, & Cassen 1984; Shu, Adams, & Lizano 1987). On the 
other hand, a powerful outflow is accelerated and colhmated via a magnetohydrodynamic mecha- 
nism from the surface of the accretion disk (Pudritz & Norman 1983; Heyvaerts & Norman 1989; 
Shu et al. 1994; Konigl & Pudritz 2000). Infalling, outflowing and rotational motions are, thus, 
simultaneously taking place in extremely young protostar envirormients, making the morphology 
and kinematics of such regions very complex. Observations of the milhmeter thermal dust con- 
tinuum and high-density molecular tracers can, in principle, reveal the phenomena taking place 
in the innermost parts of dense cores (e.g. Saito et al. 1996; Ohashi et al. 1997, 1999; Lai & 
Crutcher 2000; Hogerheijde 2001). However, the strong interaction between the outflow and the 
surrounding dense material, which has been actually observed near the driving source in a few 
cases (e.g. L1551-IRS5: Fuller et al. 1995, Ladd et al. 1995; L1157: Gueth et al. 1997; L1527: 
Ohashi et al. 1997, Hogerheijde et al. 1998, Motte & Andre 2001), complicates this portrait. In 
order to get an accurate picture of the complex envirormient of Class objects and to interpret the 
gas kinematics in these sources, high angular resolution interferometric observations of different 
tracers are definitely required. 

L1157-mm is a Class object located at 440 pc with a iboi — 11 Lq. It is associated with 
IRAS 20386+6751 and drives a spectacular outflow. The LI 157 outflow has been studied in 
details through many molecular lines, such as CO (Umemoto et al. 1992; Gueth, Guilloteau, & 
Bachiller 1996; Bachiller & Perez Gutierrez 1997; Hirano & Taniguchi 2001), SiO (Zhang et 
al. 1995, 2000; Gueth, Guilloteau, & Bachiller 1998; Bachiller et al. 2001), H2 (Hodapp 1994; 
Davis & Eisloffel 1995), NH3 (Bachiller, Martm-Pintado, & Fuente 1993; Tafalla & Bachiller 
1995), or CH3OH (Bachiller et al. 1995, 2001 ; Avery & Chiao 1996). Many other lines have been 
detected (Bachiller & Perez Gutierrez 1997; Bachiller et al. 2001), making LI 157 the prototype 
of chemically active outflows. Regarding the protostar itself, dust continuum observations have 
been carried out at 2.7 mm (Gueth et al. 1996, 1997), 1 .3 mm (Shirley et al. 2000; Motte & Andre 
2001; Chini et al. 2001; Gueth et al. 2003), 850 /xm (Shirley et al. 2000; Chini et al. 2001), and 
450 ^m (Chini et aL 2001). 

In a first attempt to understand the morphology of this young object, Gueth et al. (1997) 
carried out ^''CO (J=1^0), C^^^O (J=1^0), and 2.7 mm continuum IRAM Plateau de Bure 
interferometric observations of the region around L1157-mm, with a ~2'.'5 angular resolution. 
The continuum emission clearly shows two components: a marginally resolved compact, flat- 
tened core, perpendicular to the outflow direction, and a low-level extended emission that seems 
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to delineate the heated edges of the cavity excavated by the CO outflow. The ^■'CO emission 
also originates from the limb-brightened edges of the outflow, whereas the C^^O emission is 
more directly associated with the compact continuum source, and shows marginal evidence of 
rotation. In addition, redshifted self-absorption is present in the ^'^CO spectrum, which suggests 
infall motions. However, in order to confirm the picture described by Gueth et al. (1997), higher 
angular resolution and higher frequency continuum observations, together with different, higher 
density, tracers were necessary. In this work we present these new interferometric observations 
of file dust continuum at 2.7 and 1.3 nun, and of flie HC3N (J=12^11) and C^^O (J=2^1) 
emission around L1157-mm. These new maps reveal important details of the structure of the 
central source, and of the interaction between the outflow and the high-density core. 

2. Observations 

2. 1. Interferometric data 

Observations were carried out with the IRAM Plateau de Bure Interferometer (Guilloteau et 
al. 1992) between November 1995 and April 1996. Five different configurations of the four- 
antennas array were used. The longest baseline was 288 m (B2) and the shortest 24 m (B2, D). 
The phase center was set to the position a(J2000) = 20'^39'"06!19, (5(J2000) = 68°02'15'.'9, 
which is the L1157-mm position reported by Gueth et al. (1997). The source was observed si- 
multaneously at 109.2 GHz (HC3N J=12^11) in USB and 219.6 GHz (C^^O J=2^1) in LSB 
band. The spectra were analyzed using a correlator with one band of 10 MHz centered on the 
HC3N ( J=12^1 1) fine, one band of 20 MHz centered on the C^^O ( J=2^1) line, and two bands 
of 160 MHz for the 2.7 and 1.3 mm continuum. The spectral resolution was 0.1 1 kms^^ in both 
narrow bands. The bandpass of the receivers was calibrated by observations of the strong quasars 
3C273, 3C454.3, 0923H-392, or 1823H-568. AmpHtude and phase calibrations were achieved by 
monitoring 1823H-568, whose flux density was determined relative to 3C273 and 3C454.3. 

However, the emission of 1823H-568 turned out to be linearly polarized, and therefore pre- 
sented a sinusoidal dependence of the amplitude with time i.e. with parallactic angle. Special 
considerations had thus to be taken during the calibration process. Since the polarization of the 
two receivers (2.7 and 1 .3 mm) is orthogonal, the amplitude curve A{t) of a linearly polarized 
source has an opposite behavior on both receivers. This effect is only noticeable in tracks of 
several hours. To fix this polarization problem, we first ran a procedure to fit the calibrator am- 
plitude curve with a sinusoidal function of the parallactic angle ^I*, for each scan. The functions 
were proportional to 1h-C1 x cos^(5' + for the 2.7 mm receiver, and to 1h-C2x sin^(5' + <I>) 
for the 1.3 mm receiver, where CI and C2 give an estimate of the percentage of polarized flux, 
and $ gives an idea of the polarization angle. Then, we calibrated the amplitude curve against 
one of the well-known calibrators also used in the observations (e.g. 3C273), in order to obtain 
the real value of the flux for each scan. Thus, instead of a constant value for the flux of the am- 
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plitude calibrator, we obtained a function of the parallactic angle 5*, and this function is the one 
we used in the calibration process. 

We estimate that the uncertainties on the absolute flux density calibration are '--^20% at 
2.7 mm, and ^40% at 1.3 mm, due to the difficulties encountered to calibrate the amplitude. 
Typical rms phase noise was better than 10° at 2.7 mm, and better than 20° at 1.3 mm. The 
data were calibrated and analyzed with the GILDAS software package developed at IRAM 
and Observatoire de Grenoble. The synthesized CLEANed beam for maps made using natural 
weighting was 2'.'24 x 2'.'04 at P.A. = 33° at 2.7 mm, and 1'.'24 x 1'.'14 at P.A. = 71° at 1.3 mm. 
Unless explicitly mentioned, we subtracted the continuum from the line emission. This was per- 
formed directly in the uv plane in order to avoid non linearity effects in the deconvolution, and 
thus any ampUfication of errors induced in this process. 

2.2. Bolometer data 

Gueth et al. (2003) presented 1.3 mm bolometer measurements obtained with the MPIfR 19- 
channel bolometer array (MAMBO) on the IRAM 30-m telescope. We used this map to extract 
short-spacings information and thus complement the interferometric data set. The interferometer 
and 30-m data at 1.3 mm were merged in the uv plane using the procedure available in the 
GILDAS package. The relative weights of the interferometer and single-dish data were adjusted 
so fliat the weight density of the single-dish data matches that of the shortest interferometric 
baselines. A multiplicative scaUng factor of 0.8 was applied to the single-dish data amplitudes, 
in order to compensate for the fact that the single-dish observations were carried out at a central 
wavelength of ^1.25 mm, whereas the interferometric observations were indeed done at 1.4 mm. 
Hereafter, we will refer to the combined interferometer plus single-dish 1.3 mm continuum data 
simply as 1.3 mm continuum data. The synthesized (natural weight) CLEANed beam for the 
resulting map, 1'.'24 x 1'.'22 at P.A. = 27°, is very similar to the beam without the bolometer 
data. 

3. Continuum emission 

The millimeter continuum emission towards the core of LI 157 was mapped at 2.7 and 1.3 mm. 
Figure^shows the uniform-weight map of the 2.7 mm emission, and the natural-weight map of 
the 1.3 mm emission. Figure |2l shows the superposition of the integrated ^^CO ( J=1^0) emis- 
sion (adapted from Gueth et al. 1996), which is tracing the molecular outflow, with the 2.7 mm 
and 1.3 mm natural-weight maps (as well as with the HC3N and C^^O integrated emission). The 
continuum dust emission is resolved at both wavelengths and clearly shows two components, a 
compact source plus an extended envelope. In Table ^ we summarize the observed properties 
of these two components. The position found for the compact source LI 157-mm at both wave- 
lengths is a(J2000) = 20''39"^06?24, 5(12000) = 68°02'15'.'6. Position accuracy is about 0'.'2. 
This is in agreement with the position found by Gueth et al. (1996) through previous 115 GHz 
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Fig. 1. (a) 2.7 mm continuum emission uniform-weight map. The synthesized beam is 1'.'70 x 
1'.'33 at P. A. — 38°, and is drawn in the bottom left corner. The rms noise of the map 
is 0.33 mJybeam"^. The contour levels are —1, 1, 3, 5, 7, 9, 13, and 17 mJybeam^^. (6) 
1.3 mm continuum natural- weight map after adding the single-dish data. The synthesized beam 
is 1'.'24 X 1'.'22 at P. A. = 27°, and is drawn in the bottom left corner. The rms noise of the map 
is 1.4 mJybeam^^. The contour levels are —5, 5, 10, 15, 20, 30, 40 and 50 mJybeam^^. The 
circle represents the Plateau de Bure primary beam (50% attenuation level). 



observations. The position obtained by Gueth et al. (1997) is slightly different, but this was traced 
back to to an incorrect value for the coordinates of the phase calibrator 2013H-370. 



3.1. 2.7 mm continuum 

The 2.7 mm uniform- weight map (Fig.Q marginally resolves the compact component: the cen- 
tral source is slightly elongated perpendicular to the outflow direction, but this elongation is how- 
ever dominated by the beam. The deconvolved size we found after fitting an elliptical Gaussian 
to the uniform-weight map is ~1'.'3 x 1'.'2 at P.A. ~ 50°, which is different from the position 
angle of the uniform-weight synthesized beam (P.A. = 38°). Fitting an elliptical Gaussian di- 
rectly to the visibility data, for baselines longer than 80 m, we found a total flux of 25 mJy and a 
size of 1'.'2 ± G'.'04 x 0'.'9 ± O'.'OS at P.A. = 50°. This size, which coiTesponds to a Hnear size of 
^500x400 AU at the distance of the source, is consistent with the previous determination from 
Gueth et al. (1997). However, the flux density value that we have found is slightly lower than the 
value of 35 mJy measured by these authors. 

Extended emission at a low 1 to 3 mJybeam^^ level is also visible along the outflow axis 
at P.A. of ^-^155°. This emission has a flux density of ^31 mJy integrated over an area of 
^200 arcsec^. It is less extended and weaker than the emission detected by Gueth et al. (1997) 
with a similar synthesized beam. Adding the compact component, the total continuum flux at 
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Fig. 2. Overlay of the CO (J=l— >0) integrated emission {greyscale; adapted from Gueth et 
al. 1996) and (a) the 2.7 mm continuum natural-weight map toward LI 157 {contours). The rms 
noise of the map is 0.30 mJybeam"^. The contour levels are —1, 1, 3, 5, 7, 9, 13, 17, and 
2 1 mJy beam" ^ ; (6) the 1 .3 mm continuum natural- weight map after adding the single-dish map 
data {contours). The rms noise of the map is 1.4 mJy beam~^. The contour levels are —5, 5, 10, 
15, 20, 30, 40 and 50 mJybeam"^; (c) the HC3N (J=12— >11) integrated emission {contours). 
The rms noise of the map is 4.8 mJybeam"^ km s~^ (0.11 K km s~^). The contours are —12, 
12, 24, 36, 48, 60 and 72 mJy beam~^ km s~^ and {d) the C^^O (J=2^1) integrated emission 
{contours). The rms noise of the map is 16 mJy beam~^ km s~^ (0.28 K km s~^). The contours 
are -70, 70, 140, 210, and 280 mJybeam-^ km s'^ The HC3N and the C^^O emission have 
been integrated over the outflow velocity range 1.8 to 3.6 km s~^. The black Unes outline the 
edges of the ^^CO outflow, for reference in other figures. The circles represent the Plateau de 
Bure primary beam (50% attenuation level). 

2.7 mm is ^56 mJy. This is lower than the flux of ~90 mJy measured by Gueth et al. (1997). 
These discrepancies reflect the difficulty of absolute flux calibration in the mm domain, as well as 
that of low-level extended structures deconvolution - which may crucially depend on the actual 
uv coverage. 
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Table 1. Properties of the compact and extended continuum emission 





Integrated fluxes 


Source size'^ 




Spectral indexes 
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Component 


2.7 mm'' 
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(mJy) 


(mJy) 


(mJy) 


(") 


(AU) 


P.A. 


a 


= a - 2 


(Mo) 


COMPACT 


35 


25'' 


78" 


1.05 X 0.58' 


~450x250 


12° 


~2.1 


0.1 


~0.12 


EXTENDED 


~60 


~31 


~680 


~8x7 


-3500x3100 


~155° 


4-5 


2-3 


~1.08 


T0TAL«5 


95 


56 


760 












~1.2 



Peak coordinates of the compact core: q(J2000) = 20''39™06.24' 5(12000) = 68°02'15.6". Position accu- 
racy is0'.'2. 

(a) Integrated flux density measured by Gueth et al. (1997). 

(b) Integrated flux density measured in this paper, corrected for primary beam response. 

(c) Measured at 1.3 mm. 

(d) Integrated flux density measured by fitting an elliptical Gaussian directly to the visibility data for 
baselines > 80 m. 

(e) Integrated flux density measured by fitting an elliptical Gaussian directly to the visibility data for 
baselines > 40 m. 

(f) Size obtained by fitting an elliptical Gaussian directly to the visibility data. The uncertainty in the size 
values is ~0'.'04. 

(g) Compact + extended components. 

3.2. 1.3 mm continuum 

At 1.3 mm, the compact component has a peak intensity of 60 mJybeam^^, and is marginally 
resolved, with a size of ^I'.'l x (f.'9. By fitting an elliptical Gaussian to the visibility data for 
baselines longer than 40 m (which is consistent with using 80 m at 2.7 mm in order to select 
the same physical scales), we obtained a total flux of ~78 mJy, and a size of I'.'OS ± 0'.'04 x 
0'.'58 ± 0'.'04 (-450x250 AU) at P.A. = 12°. Note that both the size and orientation of the 
compact source are slightly different at 2.7 and 1.3 mm, which is probably related to the source 
not having a perfect Gaussian shape structure. In particular, at 1.3 mm the compact component 
is slightly elongated along the outflow cavity. 

Figures Hand 121 show that the centrally peaked compact source is surrounded by a structure 
extended roughly along the flow axis. The flux density is 5 to 15 mJybeam^^ in this component, 
well above the rms noise, ~ 1.4 mJybeam^^ (estimated over an empty area of the map). 
We measured a total flux density of ~760 mJy (corrected for primary beam attenuation). This 
is consistent with the flux of 630 mJy measured at 1.1 mm by Motte & Andre (2001) with 
MAMBO on the IRAM 30-m telescope. The circumstellar dust emission around L1157-mm is 
clearly dominated by the extended component. This envelope does not have spherical symmetry, 
but it is spread over a region of —8" (~3500 AU) along the outflow direction. 

The 1.3 mm map also shows an emission extended in the direction perpendicular to the 
outflow, on a ^-^20" scale. This structure is obviously much too large to be a circumstellar disk. It 
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Fig. 3. Velocity channel maps of the HC3N (J=12^11) emission. The wlsr of L1157-mm is 
2.6 km s^^. The central outflow velocity of each velocity interval is indicated at the upper left 
corner of each panel. The l-cr noise in 1 channel is 7.5 mJybeam^^. Contour levels are —15, 
15 to 75 by step of 15 mJybeam^^. The conversion factor from Jybeam^^ to K is 22.5. The 
synthesized beam is 2'.'24 x 2'.'04 at P. A. = 33°, and is drawn in the bottom left corner of the 
bottom left panel. The dash lines outline the edges of the outflow (see Fig.O. 



might be due to a flattened envelope remnant of the molecular cloud in which the protostar was 
formed. This extended material has also been detected, but at a much larger spatial scale, through 
850 fim and 1.3 mm single-dish observations (Shirley et al. 2000; Chini et al. 2001; Gueth et al. 
2003). 

4. Molecular line emission 

4.1. HC^Nmaps 

The HC3N (J=12^11) transition has a high critical density (8 x 10^ cm^'', Chung, Osamu, & 
Masaki 1991), which makes it an even better tracer of high density regions than e.g. CS ( J=1^0) 
and CS ( J=2^ 1). Figure|3]shows the velocity channel maps for the HC3N ( J=12^ 1 1) emission 
toward the core of LI 157, around the systemic velocity channel (~ 2.6 km s^^). The gas emis- 
sion around LI 157-mm is not compact, but extended and elongated along the outflow edges. The 
outflow is blueshifted toward the south and redshifted toward the north and is almost in the plane 
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Fig. 4. Velocity channel maps of the C^'^O (J=2^1) emission. The wlsr of L1157-mm is 
2.6 km s^^. The central outflow velocity of each velocity interval is indicated at the upper left 
corner of each panel. The l-tr noise in 1 channel is 30 mJybeam^^. Contour levels are —50, 
50 to 200 by step of 50 mJybeam^^. The conversion factor from Jy beam^^ to K is 17.9. The 
synthesized beam is 1'.'24 x 1'.'14 at P. A. = 71°, and is drawn in the bottom left corner of the 
bottom left panel. The dash lines outline the edges of the outflow (see Fig.|2ji. 

of the sky (Gueth et al. 1996). As a consequence, there is blueshifted and redshifted HC3N emis- 
sion in both lobes. The integrated HC3N emission (Fig.|2j has a flux density '^670 mJykms^^. 

At a distance of '^30" from the central source, the southern lobe reveals an elongated fila- 
ment at velocity 2.6 and 2.8 kms^^ (Fig.|3}. It is outlining very clearly the eastern flank of the 
outflow CO lobe. This emission is already significantly affected by the primary beam attenua- 
tion. Interestingly, this feature coincides with a SiO shock (S3 shock from Gueth et al. 1998; 
Zhang et al. 2000), traced also by H2 (Davis & Eisloffel 1995), NH3 (Tafalla & Bachiller 1995), 
H2CO, CS, CH3OH, and SO (Bachiller et al. 2001). Gueth et al. (1996, 1998) and Zhang et 
al. (2000) describe a precessing episodic jet scenario that can explain the different orientations 
and velocities of the cavities of the outflow. According to this scenario, the jet has precessed 
toward the SE-NW direction, and the region discussed here would thus correspond to an event 
of enhanced ejection, which is impacting the flank of the existing CO cavity. The resulting shock 
explains the increased abundances of some species at that location, including HC3N (Bachiller 
& Perez-Gutierrez 1997). The velocity of the HC3N gas is significantly lower than that of the CO 
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emission at the same location (e.g. Gueth et al. 1996). This suggests that high velocity material 
has a lower density, and would thus not be seen in HC3N, although we cannot rule out an opacity 
effect. No HC3N emission has been detected inside the cavity, but an extended emission could 
possibly have been filtered out by the interferometer. 

4.2. C^^Omaps 

Figure 13 presents the velocity channel maps for the C^'^O (J=2^1) emission. The integrated 
emission (Fig.|2j appears more compact than the HC3N emission: it is dominated by the central 
compact component, and is thus associated to the innermost part of the protostellar envelope. 
The total integrated intensity is ~3 Jy kms^^. The C^^O emission is also tracing gas associated 
to the outflow, as can be seen at the systemic velocity channel map (ulsr= 2.6 kms^^, see 
Fig-E}- A. striking feature visible in Fig.|4]is that the C^^O emission at 2.8 kms^^ is weaker 
than in the previous and following velocity channels. This effect is clearly seen as a narrow deep 
self-absorption in the spectrum taken at the position of LI 157-mm (see below, section l6.2.1l and 
Fig.EJ. 

5. Analysis 

5. 1. Spectral index map 

Figure |5]presents a map of the spectral index a, where 5*^ oc ly", computed using the 2.7 and 
1.3 mm continuum interferometric images restored with the same resolution, 2'.'2 (or ^900 AU). 
The spectral index varies from a value ~2.1, at the position of the compact source to a value 
of ^^4-5 at larger distance'. Instrumental effects such as missing flux and/or a wrong flux scale 
would affect the absolute level of the spectral index map, but can hardly introduce a spatial 
variation. Line contamination can be excluded, since it does not affect the interferometer data 
and was found to be absent in the bolometer measurement (Gueth et al. 2003). The observed 
variation accross the L 1 157-mm region may thus be indicative of physical differences between 
the properties of the material toward the embedded source and those of the more extended region. 
This effect was already observed in other young embedded sources (Visser et al. 1998; Johnstone 
& Bally 1999; Smith et al. 2000; Beltran et al. 2002). 

A low value of a, as observed at the central position, can be due to optically thick dust 
emission. In such a case, the dust temperature Td is equal to the brightness temperature (in 
the Rayleigh- Jeans approximation), which can be derived from the observed flux density: 



' Gueth et al. (2003) have analyzed the spectral index around LI 157-nim from bolometric observations at 
1.2 and 0.85 mm. They obtained an intermediate value of a ~ 3.5 in the central 20" beam which includes 
all our interferometric map. 
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where fl is the sohd angle subtended by the source. An upper limit on Td can be computed 
when assuming that the whole flux of the compact core (78 mJy, see Table 1) is coming from 
the optically thick contribution. Assuming that the size of the optically thick emitting region is 
around 150 AU, we derive a value of < 10 K for the dust temperature (note that the larger the 
optically thick core, the lower Tj). This value is very low and suggests that a significant fraction 
of the compact core has an optically thin emission. 

Assuming that the dust emission at mm wavelengths is optically thin and in the Rayleigh- 
Jeans regime, a is related to the power-law index /3 of the dust emissivity oc v^, through 
a = /? + 2. For the more extended material, /3 would thus be larger than 2, the typical value for 
interstellar dust grains (see e.g. Draine & Lee 1984), while the value of /3 for the compact source 
would be -^0.1. Interestingly, this decrease of (3 could be due to larger grains (e.g. Mannings & 
Emerson 1994), hence suggesting grain growth in the denser inner regions. Grain shape evolution 
or chemical evolution are also possible (see, e.g. Ossenkopf & Herming 1994; Pollack et al. 
1994). 




20^39"'07'0 06^0 05^0 

a(2000) 



Fig. 5. Spectral index map from the 2.7 and 1.3 nmi continuum images, restored to the same 
resolution of 2'.'2. Contour levels are plotted from 2 to 5 by intervals of 0.5. Greyscale is Unear in 
the range of 2-5, and shows the regions with higher spectral index darker. The cross marks the 
position of LI 157-min. 
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5.2. Envelope mass - from continuum emission 

Assuming optically thin dust thermal emission, the mass of the emitting material can be derived 
from the 1.3 mm continuum flux density Si, by: 

Using a dust temperature of 40 K (following Gueth et al. 1997), and a dust mass opacity of 
0.01 cm^g^^, as recommended by Henning, Michel, & Stognienko (1995) for a gas-to-dust 
ratio of 100, we obtain a mass of the central compact source of ^-^0.12 Mq, while the mass of the 
envelope is ^^1.08 Mq. Thus, the total (compact + extended) circumstellar mass is ^^1.2 Mq. 
This is a factor 3 below the value of ^3Mq derived by Gueth et al. (1997) at 2.7 mm. The 
absolute flux density of the 1.3 mm data may be significantly underestimated. These calculations 
are also dependent on the dust opacity law used, the mass being lower in case of higher dust 
opacity. In any case, the large ratio between the mass of the extended component and that of the 
compact one indicates that there is an important reservoir of material that can be incorporated 
into the central protostar. This is in agreement with the conceptual definition of a Class object. 



5.3. Envelope mass - from gas emission 

The gas mass of the protostellar condensation can be better estimated from the C^**0 data, since 
the HC3N emission is partly associated to the outflow. Following the derivation of Scoville et al. 
(1986), and assuming optically thin emission, the C^^O beam averaged column density is given 
by 

N = 1.15 IQi^ + J Tb dv (3) 

where Tex is the excitation temperature, and / Tsdv is the integrated brightness temperature 
of the C^*0 ( J=2— >1) emission. Assuming Tex = 40 K, a C^^O abundance relative to molec- 
ular hydrogen of 1.4 x 10^^, and a mean molecular weight of 2.6 mn, we obtain a mass of 
^9x10^'^ Mq in the inner ^2" x 2" area (roughly corresponding to the central source). It 
should be mentioned that the excitation temperature value adopted does not affect significantly 
the estimated molecular mass; the mass is ^7 x 10~^ Mq for Tox = 20 K, and ^ 1 1 x 10"-^ Mq 
for Tox = 50 K. Our derived value is certainly a lower limit since the C^^O is affected by self- 
absorption (see 5 16.2. It . and hence partly optically thick. Indeed, Gueth et al. (1997) derived a 
mass of 4.7 x 10~^ Mq fromC^^O (J=1^0). Comparing the molecular mass, a few 10^^ Mq, 
to that derived from the continuum data, M ~ 0.12 AIq, it suggests that CO in LI 157-mm could 
be depleted by a factor around ^3-10, although this result has to be taken with caution in case 
the gas is optically thick and/or clumpy. Using the same assumptions as above, the mass of the 
whole C^^O envelope (~8" x 8" area) is ~4xlO^^ Mq, again significantly below the mass 
derived from the continuum emission (1.08 Mq). However, the C^^O emission is less extended 
that the continuum emission, and the comparison between the two masses may therefore not be 
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Fig. 6. HC3N (J=12^11) (left) and C^^O (J=2^1) (ng/jf) spectra obtained at the position of 
L1157-mm. The continuum emission has not been subtracted (the continuum level is indicated 
by the dashed horizontal line). The dashed vertical line indicates the ulsr of 2.6 km s^^. The l-a 
noise in 1 channel is 0. 17 K for HC3N, and 0.45 K for C^^O. The conversion factor of Jy beam~^ 
to K is 22.5 for HC3N, and 17.9 for C^^O. 

relevant. It is also possible that we miss some C^^O ( J=2-^l) emission because of the lack of 
short spacings in our measurement. 



6. Discussion 

6. 1. Interaction between the outflow and the envelope 

Both our continuum and lines observations reveal structures or filaments that seem to be associ- 
ated with the edges of the L 1 157 outflow. 

The 1.3 mm continuum map (Fig.|3 reveals a low-level extended emission whose conical 
shape matches that of the outflow lobes. Gueth et al. (1997) reported a similar feature from 
2.7 mm observations: the dust emission was observed and modelled as a "cross-like" structure 
coincident with the outflow edges. Note however that these observations indicated opening an- 
gles slightly smaller than those shown by our maps. The morphological association between the 
outflow edges and dense protostellar material in L 1 157 is even better shown by the the HC3N 
observations. The channel maps presented in Fig.|3l(see the 2.6 to 3.2 km s^^ channels) reveal 
extensions that outline quite accurately the beginning of the outflow conical lobes. The south- 
eastern arm of this "cross-like" pattern is missing, as in the 2.7 mm continuum from Gueth et 
al. (1997), which suggests different temperature and/or density than in the other arms. The C^^O 
emission, while much more compact than the HC3N emission, exhibits marginal indications of a 
similar morphology. 

Following Gueth et al. (1997), we suggest that these morphological coincidences indicate 
that a strong interaction between the outflow and the protostellar envelope is taking place in 
L 1157: the outflow sweeps up envelope material, thereby creating heated/compressed regions 
at its edges, which, observationaly, mimic a "cross-like" structure (see Gueth et al. 1997, their 
Fig. 10). Circumstellar dust and gas material around young stellar objects outlining the edges of 
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the molecular outflow have also been mapped in other sources (e.g. B5-IRS1: Langer, Velusamy, 
& Xie 1996; L1527: Ohashi et al. 1997, Hogerheijde et al. 1998, Motte & Andre 2001; L1551- 
IRS5: Fufler et al. 1995, Ladd et al. 1995). 

Interestingly, the HC3N emission also revels two peaks located on both sides of the central 
protostar. They are clearly visible in the integrated emission map (Fig. ^ as well as in the in- 
dividual channel maps (Fig. |3. Their positions coincide with the apparent base of the outflow 
conical cavities as traced by the ^^CO (see Fig.|2}, suggesting that these emission peaks could be 
produced when the dense gas is shock-heated and compressed by the flow. Hence, we conclude 
that, at the resolution of our observations, the structure of the dust-ngas protostellar envelope is 
strongly affected by the presence of the outflow. 

6.2. Kinematics of ttie Gas 
6.2.1. C^^O Self-absorption 

Figure|6|presents the HC3N (J=12^1 1) and C^*0 ( J=2-^l) spectra obtained at the position of 
LI 157-mm. As already mentioned, the C^^O emission shows a narrow (At; ~ 0.2 km s^^) deep 
absorption feature in the spectra, at a velocity of ^--^2.8 km s~^. A similar feature of ^0.3 km s^^ 
was also reported at the same velocity in the interferometric ^■^CO (J=1^0) map and spectra 
by Gueth et al. (1997). This suggests that the C^^O (J=2^1) line is self-absorbed. However, 
missing short-spacings information could cause a fake absorption, but spatially extended (see 
Gueth et al. 1997, their Fig.7). In Fig.Qwe plot the C^^O ( J=2^1) spectra ai'ound LI 157-mm, 
by offsets of 1'.'2, roughly the size of our synthesized beam. This figure shows that the absorption 
feature is only detected at a few positions, close to the compact source and not at larger distance 
of LI 157-mm. Hence, the self-absorption, confined in a small region (< 5" or ^2000 AU) 
around LI 157-mm, appears real. Note that Fig.0also shows line broadening localized toward 
the center 

A redshifted self-absorption is naturally explained in a scenario where an optically thick 
line is observed in an infalling core with temperature decreasing outward (e.g. Leung & Brown 
1977). Fig.|6lshows that the brightness temperature (in T,„b scale) of the self-absorbed feature 
is ^2.2 K. Assuming optically thick emission, it corresponds to an excitation temperature of 
Tox — 6.3 K, in excellent agreement with the excitation temperature found by Gueth et al. (1997) 
from the ^'^CO spectra. This low Tex indicates that the self-absorption features are most probably 
sub-thermally excited, and should thus be produced in a medium whose density is low enough 
to prevent thermalization of the lines. No self-absorption was detected in the C^^O (J=1^0) 
interferometric maps obtained by Gueth et al. (1997). This suggests that the (J=l-^0) line is 
optically thin while the (J=2^1) is significantly thicker. For Tox ~ 10 K, the ratio of optical 
depths between the two transitions is t(J=2^1)~ 1.9 t( J=1— +0). Note also that the C^^O 
(J=1^0) spectra from Gueth et al. (1997) has a narrow profile, making it difficult to detect a 
possible asymetry. 
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Fig. 7. C^^O (J=2^1) spectra obtained around the L1157-mm central region. The continuum 
has not been subtracted from the line emission. Offsets from the central position are indicated (in 
") at the upper left corner of each plot. The dashed vertical line indicates the wlsr of 2.6 km s^^. 



6.2.2. Other absorption features? 

In addition to the self-absorption feature discussed in the previous paragraph, Fig.|6lreveals that 
both the HC3N and C^^O spectra at the position of L 1 157-mm have weaker emission in the 1- 
2 km s^^ (blueshifted) and the 4-5 km s^^ (redhifted) velocity intervals. This suggests that the 
actual line profile may be triangular-shaped, with two strong absorption features. In both lines, 
each of the "absorbed" channels is at a few cr below the un-affected emission, but the whole 
absorption feature is significant because it extends over 3 to 6 adjacent channels. The redshifted 
feature in the C^^O even shows an aborption of the continuum emission. 

To explain such absorptions in the spectra, one has to invoke the presence of cold absorbing 
gas in front of a warm emitting region. Concerning the blueshifted feature, it could arise from 
self-absorption in the blueshifted southern outflow lobe (see Fig.lSjl: cold gas at the cavity edge 
is located in front of warmer material present along the flow axis; the latter gas is expected to 
have larger velocities, but since it is closer to the plane of the sky, projection effects may produce 
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Fig. 8. Geometry of the disk/outflow system in L 1157 and possible origin of the blueshifted 
absorption feature. 

the same line-of-sight velocity for both components. As for the redshifted absorption feature, a 
similar mechanism in the redshifted, northern lobe cannot be proposed because the cold material 
would be located behind the warm gas (Fig. |8|l. Moreover, any phenomenum in the redshifted 
outflow lobe could not explain the absorption of the continuum, since this lobe is behind the 
central source (which is at the origin of the compact continuum emission). Instead, one could 
invoke absorption by the infalling envelope, whose redshifted part is in front of the continuum 
source. 

This latter scenario can drive interesting constrains on the size of the continuum emitting 
region. The velocity at which the redshifted continuum absorption is observed is higher than the 
self-absorption feature: this means than the absorbing gas is located closer to the center, and thus 
should be hotter The excitation temperature should thus be higher than 6.3 K. On another hand, 
the excitation temperature must be significantly smaller than the continuum brightness temper- 
ature (1 K), in order to create an absorption. These constrains can be reconciled if we assume 
that the continuum emission is dominated by a strongly peaked central component, which is af- 
fected by beam dillution effects. Let's assume that the C^'^O absorbing material has an excitation 
temperature of 10 K. Judging from the absorption feature (Fig.|6|l, the continuum true brightness 
temperature must be typically ^5 times larger, i.e. ^^50 K. This can be explained by a beam 
dilution of '^50 in area, hence ^1 in size. This in turn would imply a continuum source size of 
^75 AU (beamsize/7~ 0'.'17). Interestingly, this value is of the order of the expected typical size 
of a circumstellar disk in a Class object (e.g. Pudritz et al. 1996). 
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Clearly, these absorption features have to be confirmed and would deserve further, more 
detailed, observations in order to derive a clear understanding of the complex kinematics of the 
LI 157 protostellar envelope. 

6.2.3. Rotation ? 

It should be mentioned that we did not find any evidence of velocity gradient for the 
C^^O ( J=2— >1) transition (see Fig.|4}, indicative of rotation, as suggested by Gueth et al. (1997) 
fortheCi*O(J=1^0)line. 

7. Summary 

We have studied with the IRAM Plateau de Bure millimeter interferometer the dust and gas 
emission toward the core of the Class object L1157-mm. Our main conclusions are: 

- The continuum emission at 2.7 and 1.3 mm show a compact component which is resolved 
at 1.3mm into a flattened compact core of ^450x250 AU, and mass ^^0.12 Mq and is 
surrounded by an extended envelope of ^3000 AU, and mass ~ 1 . 1 Mq . 

- The milUmeter spectral index varies across the region. A lower value a 2.1 is found to- 
ward the position of the compact protostar, while the spectral index is a > 4 for the extended 
surrounding material. These values imply a dust emissivity index ~ 0. 1 and /3 > 2, respec- 
tively. This variation could possibly indicate grain growth toward the compact component. 

- A strong interaction between the molecular outflow and the close protostellar environment is 
taking place and affects the structure of the innermost parts of the envelope. This is shown by 
the spatial correlation between the molecular outflow and the dust continuum (1.3 mm map) 
and the HC3N emission: both maps show structures associated to the edges of the outflow 
lobes, as traced by the ^^CO emission. 

- HC3N emission is also detected at a distance of ^13000 AU from the central object, at the 
position of a shock associated to the impact of the precessing jet against the walls of the main 
cavity of the southern lobe. 

- Evidence of infall has been detected through the C^^O ( J=2^1) observations, in agreement 
with the indications of infall detected through ^'^CO by Gueth et al. (1997). 

- We did not detect any velocity gradient indicative of a rotating circumstellar disk. 

Finally, when observing Class environment such as that of L 1 157-mm, multiple velocity 
components and temperature gradients along the Une of sight do confuse the interpretation of 
spectral line emission. In particular, they make it difficult to disentangle the various components 
(disk, envelope, outflow) and determine their velocity distributions. Our results show that even 
the highest density gas tracers such as HC3N are affected by the ejection phenemena, tracing 
also ambient material compressed by the outflow. When deriving results from low angular res- 
olution observations of Class objects one should be aware of these problems. High angular 
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resolution observations of different dust and gas tracers are absolutely necessary to get a detailed 
and accurate picture of such young objects environment. 
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